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ABSTRACT 

An analytical formula is developed to accurately represent the photoabsorption cross section of O i for all energies 
of interest in X-ray spectral modeling. In the vicinity of the K edge, a Rydberg series expression is used to fit 
/x’- matrix results, including important orbital relaxation effects, that accurately predict the absorption oscillator 
strengths below threshold and merge consistently and continuously to the above-threshold cross section. Further, 
minor adjustments are made to the threshold energies in order to reliably align the atomic Rydberg resonances 
after consideration of both experimental and observed line positions. At energies far below or above the K-edge 
region, the formulation is based on both outer- and inner-shell direct photoionization, including significant shake-up 
and shake-off processes that result in photoionization-excitation and double-photoionization contributions to the 
total cross section. The ultimate purpose for developing a definitive model for oxygen absorption is to resolve 
standing discrepancies between the astronomically observed and laboratory-measured line positions, and between 
the inferred atomic and molecular oxygen abundances in the interstellar medium from xstar and spex spectral 
models. 
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1. INTRODUCTION 

Atomic photoionization, an important astrophysical process, 
has been studied for more than a century — since the seminal 
understanding of its energetics by Einstein (1905) and the 

first calculations of quantum mechanical cross sections (Bates 

1 939). Over the years, a plethora of experimental and theoretical 

investigations have managed an excellent grasp of its physics 

(Fano & Cooper 1 968 ; Starace 1982), together with a remarkable 
quantitative description of the valence-shell photoionization of 

atoms and atomic ions (Opacity Project Team 1995, 1997). 
However, the quantitative model of inner-shell photoabsorption 
is less sound because of a variety of relaxation processes, namely 
Auger and X-ray emission, that must be taken into account to 

achieve acceptable accuracy, especially in the near-threshold 
region. 

Inner-shell photoabsorption of metals with nuclear charge 

7 ^ Z ^ 28 is directly accessible to modern X-ray obser- 
vatories, such as Chandra and XMM-Newton, and, hence, is 

of much interest in astronomy. Particularly prominent in the 
photoabsorption of the interstellar medium (ISM) are the K- 

shell features (lines and edges) of atomic oxygen, which are the 

most abundant metal and is critically important in the energetic 

and chemical evolution of the universe (Stasinska et al. 2012). 
At present, though, the unsatisfactory quantitative understand- 
ing of oxygen inner- shell photoabsorption is such that there 
exist various sets of cross sections, each one leading to different 
conclusions regarding the ionization and atomic-to-molecular 
fractions in the ISM along various Galactic lines of sight. 


The first inner-shell photoabsorption cross sections of oxygen 
reported (Henke et al. 1993; Verner et al. 1993; Verner & 
Yakovlev 1995) were simple step-function fits to low-resolution 
solid-state data (Henke et al. 1993) or to theoretical calculations 
by Reilman & Manson (1979) using a central potential method. 
These results depict cross sections across inner-shell thresholds 
with unphysical, discontinuous edges where even the threshold 
energies are poorly determined. These cross sections were used 
by Schulz et al. (2002) in an early analysis of ISM absorption 
near the oxygen K edge in Chandra X-ray binary-star spectra. 

A later theoretical cross section, which took into account 
resonance effects, was computed by McLaughlin & Kirby 
(1998) using the /x’-matrix technique. However, this calculation 
failed to include the effects of orbital relaxation and spectator 
Auger damping, which causes blending of the resonances 
converging to the inner-shell thresholds, thus smearing the 
otherwise sharp K-shell edge. Both Paerels et al. (2001) and 
Takei et al. (2002) used this cross section to analyze the ISM 
K-shell absorption of oxygen in the Chandra spectra toward 
X06 14+091 and Cyg X-2, respectively. All these studies found, 
after fitting the O i K« line and edge, residual narrow absorption 
at Rs23.36 A and a broad edge feature at «s22.9 A. In all cases, 
the residual absorption was most likely attributed to oxygen 
compounds, although the narrow absorption feature could also 
be due to O ii. 

By contrast, de Vries et al. (2003), using the XMM-Newton 
Reflection Grating Spectrometer (RGS), found that the ISM 
oxygen K-shell edge observed in X-ray binaries and extragalac- 
tic sources was well described by the R-matrix cross section of 
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McLaughlin & Kirby (1998). A second A-matrix calculation of 
the O i cross section was reported by Gorczyca & McLaughlin 
(2000), with a full account of relaxation and Auger damping, 
in fairly good agreement with the laboratory measurements of 
Stolte et al. (1997). As a result, Juett et al. (2004) analyzed the 
Chandra spectra available at the time using the cross sections by 
both McLaughlin & Kirby (1998) and Gorczyca & McLaughlin 
(2000). They point out that the A-matrix cross sections and those 
by Verner et al. (1993) and Verner & Yakovlev (1995) agreed 
to within ^5% well above threshold, but with significant differ- 
ences in the threshold region. These differences were such that, 
when using the more recent cross section to fit the spectrum, the 
previously found broad threshold residuals disappeared. Thus, it 
was concluded that the narrow absorption feature in the spectra, 
after subtracting the O i contribution, was due to trace amounts 
of ionized oxygen rather than to molecular compounds. Juett 
et al. (2004) also found that the discrepancies between the var- 
ious calculations and experiments regarding the wavelengths 
of the Oi K« line and the K-shell threshold were consider- 
ably greater than the resolution of the astronomical spectra. 
Therefore, this energy dispersion may be the main source of 
uncertainty left in the atomic cross sections. 

Since then, Garcia et al. (2005) have reported A -matrix 
calculations for the whole oxygen isonuclear sequence, which, 
in the case of O i, agree to within »H0% with the near-threshold 
cross section of Gorczyca & McLaughlin (2000). Subsequently, 
Garcia et al. (2011) used these datasets to reanalyze seven 
XMM-Newton observations of the X-ray binary Sco X-l, and, 
by adjusting the absolute wavelength scale of the theoretical 
cross sections, found that the spectra were well fitted by Oi 
absorption alone with no conclusive evidence of contributions 
from any other source. On the other hand, a thorough study of 
XMM-Newton spectra toward the low-mass binary GS 1 826-238 
by Pinto et al. (2010) showed that the ISM was composed of 
a mixture of multi-phase gas, dust, and molecules; in the case 
of oxygen, its abundance was found to be 20-30% higher than 
protosolar, and at least 10% of its column density was in the 
form of molecules and dust grains. These findings have been 
confirmed by Pinto et al. (2013) in a more extensive survey 
of nine low-mass X-ray binaries where 15-25% of the total 
amount of oxygen was found to be condensed in dust. Moreover, 
in a recent examination of several Chandra spectra toward the 
low-mass binary XTE J 1 8 17 — 330, Gatuzz et al. (2013) were 
forced to shift the photoionization cross sections of Garcia et al. 
(2005) to fit the absorption lines from both the high- (Ovi, 
Ovn) and low-excitation (Oi, On, Om) plasma components. 
The discrepancies pointed out by Juett et al. (2004) regarding 
the observed and measured positions for the Kor and K/J lines in 
O i still stood. Gatuzz et al. (2013) report an oxygen abundance 
close to solar, which, in essence, dissents from the conclusions 
reached by Pinto et al. (2010, 2013). 

Laboratory measurements of the Oi K vacancy states 
(Caldwell et al. 1994; Krause 1994; Menzel et al. 1996; 
Stolte et al. 1997; McLaughlin et al. 2013), specifically the 
ls2s 2 2p 5 V 3 " resonance, also show a bothersome scatter. This 
issue will be addressed more fully in Section 2.3.4 since this 
uncertainty, and that of the observations of interstellar oxygen 
X-ray spectra, are at the heart of the remaining issue of absolute 
energy normalization. The recent experiment of McLaughlin 
et al. (2013) is similar to the earlier study by Stolte et al. (1997), 
but the entire resonance region is now covered in one contin- 
uous scan in photon energy rather than the previous piecemeal 
scans for the lower n — 2 member and the higher 3 ^ n -> oo 


Rydberg series; however, nearly identical energy positions are 
reported. The new theoretical results, on the other hand, are 
obtained from a 910-level A’-matrix calculation that, upon close 
inspection, are essentially equivalent to those in Gorczyca & 
McLaughlin (2000). Thus, there is nothing substantively new 
learned from this study; it reconfirms the same resonance en- 
ergy positions. This issue will be addressed more fully in 
Section 2.3.4. 

The ultimate purpose of the present study is to arrive at 
a consensus for the best description of the photoionization 
cross section for neutral atomic oxygen. To this end, it is ad- 
vantageous to create a single photoabsorption model that is 
transparent to all atomic data users. This is most easily accom- 
plished by formulating an analytical expression that includes all 
the desired essential features: accurate background cross sec- 
tions (the “shoulders”), line positions, widths, and oscillator 
strengths. We accomplish this by appealing to a combina- 
tion of A -matrix computations, laboratory measurements, tabu- 
lated solid-state absorption data, independent-particle (IP) data, 
multi-configuration atomic structure calculations, and astro- 
nomically observed X-ray lines. Given a consistent photoab- 
sorption cross section, we then plan to use the same atomic 
description in two different X-ray spectral modeling codes to 
render the ISM oxygen K features. 

2. ANALYTICAL MODEL CROSS SECTION 

We attempt to develop an analytical expression for the 
most reliable photoabsorption cross section possible. We begin 
with new A- matrix calculations, which are slight improvements 
over earlier results (Gorczyca & McLaughlin 2000) that were 
benchmarked favorably with experiment (Stolte et al. 1997). 
These results, together with experimental measurements, allow 
for the most accurate representation of the strong I .v — > np 
resonances below the K edge. At energies far below the K-edge 
region, where only outer- shell photoionization occurs, we use 
the formula of Verner et al. (1996), which is a simple fit to the 
IP results of Reilman & Manson (1979), and is found to be in 
excellent agreement with the present A-matrix results. At higher 
energies, a fit to the tabulated data of Henke et al. (1993) is used; 
these data are assessed to be the most accurate since important 
shake-up and shake-off processes are also accounted for (see 
Section 2.3.1). The resulting total photoabsorption cross section 
as a function of photon energy E = hv is thus partitioned as 

a PA (£) = a 2s , 2p (E) + a[ f(E) + < ect (£). (1) 

To illustrate this demarcation, we plot several data sets of the 
photoabsorption cross section in Figure 1, where the three 
energy regions are all depicted. For energies below ~520 eV, 
the cross section consists solely of the cr 2s , 2 p(E) outer-shell 
photoionization contribution, whereas, just below the K-edge 
region, the strong Is —*■ np resonance absorption profiles 
dominate. At higher energies, the cross section is essentially 
due to direct I ,v — > ep photoionization and accompanying 
photoionization-excitation and double photoionization (which 
the present Is -photoionization A-matrix calculations do not 
include as discussed in Section 2.3.1). These three regions, 
and the precise theoretical modeling of each, are henceforth 
outlined. 

2.1. Outer-shell Photoionization 

For the outer-shell photoionization cross section, we have 
verified that the fit of Verner et al. (1996) is reliable, and their 
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Figure 1. Broad view depiction of the Oi photoabsorption cross section, 
indicating where the outer-shell ( (*2s,2p X inner-shell (o r 5 rect ), and resonance 
(er[® s ) contributions are most important. Shown are the present /^-matrix (red 
curve), analytic formula (green curve), IP fit of Vemer et al. (1996; cyan curve), 
and Henke et al. (1993) data (blue squares). The jR-matrix results account for 
resonances but are missing two-electron contributions at higher energies (see 
text). The fit incorporates all the correct physics. 

(A color version of this figure is available in the online journal.) 


analytical formula 


^Is.lpiE) = CT 0 [O - if + )4] y (P U)/2 (l + yfyJVa) P (2) 

is therefore adopted, where x = E/ Eq — yo, y = Vx 2 + y and 
the fitting parameters ero, Eq, y a , p, y w , yo, and y\ are listed in 
Table 1 . The fit in Figure 1 is in close agreement with the present 
R-matrix results and with the tabulated data of Henke et al. 
(1993) at energies slightly above the 2s ionization threshold up 
through the region just below the I .v — »• np resonances. At lower 
energies, important channel-coupling effects and prominent 
outer-shell resonance structure are found in the R-matrix results 
(and modeled somewhat more crudely in Henke et al. 1993). 
However, since we are not concerned with such low energies, 
the present fit is sufficient. Just below the resonance region at 
about 520 eV, there is coupling between the (open) 2 s~ l ep and 
(closed) \s~ l np channels that gives rise to a slight dip in the 2s 
R-matrix cross section (not observable on the scale of Figure 1) 
due to a transfer of oscillator strength (Dias et al. 1997; Hansen 
et al. 1999). This small dip over such a narrow energy region is 
ignored in the present final model. 

2.2. Inner-shell (High-energy) Photoionization 

We find that the data of Henke et al. (1993), above the 1 .v 
threshold, can be accurately fitted with the expression 


a 


Henke 

Is 


(E) = (7 th 




where the fit threshold position is chosen to be = 
544.544 eV, giving a threshold cross section of <7* = 1.07 Mb 
and parameters ot\ = —0.7227 and a 2 = 0.2153, which are 
needed for further fitting (see Table 1). We choose this func- 
tional form in our effort to derive an expression that matches 
continuously from below each of the two main On thresholds, 
and, therefore, only the two parameters 04 and are needed to 
get the correct shape of the above-threshold cross section. 


Table 1 

Summary of Fitting Parameters 


Cross Section 

Parameters 

<*2s,2p 

(T 0 = 1745.0, Eg = 1.24, y a = 3.784, p = 17.64 
y w = 0.07589, y 0 = 8.698, yi = 0.1271 

°i? 

fo = 0-132, /^ = |/o, = |/o 

1v2.v 2 2/) 4 ( 4 P)np Series (i s = 1) 

E ^ = 544.54 eV, T 1 = 0.1348 eV 

= 1.11, p.\ =0.77, =0.75 (n ^4) 

/o ,2 = 0-867 f 0 \ 3 = 0.93 /g^, f 0 ‘„ = fi" ( n > 4) 
ls2s 2 2p 4 ( 2 P)np Series (i s = 2) 

El = 549.32 eV, F 2 = 0.1235 eV 
p\ = 0.84, p. 2 n = 0.80 ( n > 4) 

/o 2 ,3 = 1.02/o 2 00 ,/o 2 „ = /o 2 oo(»>4) 

-direct 

a u 

= -0.7227,0-2 = 0.2153 

Note. The fitting data for each of the two series i s = 1 and i s = 2 follow the 


respective designations (bold faced). 

2.3. Resonance Region 

Our strategy in the Is -> np resonance region is first to 
formulate an analytical fit (see Section 2.3.2) to the results from 
new A- matrix calculations (detailed in Section 2.3.1). The fit 
parameters are then adjusted slightly to match the experimental 
resonance positions and oscillator strengths, highlighting both 
the best assessment of the absolute energy scales as determined 
by line observations and the smooth consistent merging to the 
above-threshold cross section. 

2.3.1. R-matrix Calculations 

The present R-matrix approach is based closely on the earlier 
work of Gorczyca & McLaughlin (2000), the essential differ- 
ence being in the present removal of pseudo-resonances. Further 
improvements include, first, a new computation of the specta- 
tor Auger widths using a resonance time-delay matrix analy- 
sis (Smith 1960) within an independent A’-matrix calculation 
for the e~— On scattering. A second improvement involves a 
smooth turn-off of the spectator Auger damping E -> E + iT/2 
(Gorczyca & Robicheaux 1999; Gorczyca & McLaughlin 2000) 
as the effective quantum number tends to the orbital angular mo- 
mentum (v | /) to avoid the discontinuity previously seen when 
the quantum defect approach is abruptly turned off. Here, we 
allow the width to vanish continuously (T — > 0) in this limit be- 
fore the quantum defect channel is omitted (Gorczyca & Badnell 
2000 ). 

As in earlier work (Gorczyca & McLaughlin 2000), we also 
emphasize the importance of accounting for orbital relaxation 
following inner-shell photoionization: the 2 s r and 2 p r “relaxed” 
orbitals in the final l.v2.v 2 2p ( 4 Oil vacancy state differ signif- 
icantly from those in the initial Oi ls 2 2s g 2p g ground state 
because of the doubling of the effective charge seen in the O n 
state. As a result, the computation of the direct cross section in- 
volves an overlap amplitude factor proportional to the (2s g \2s r ) 
and (2p g \2p r ) orbital overlap integrals; it is, therefore, impera- 
tive to account for this orbital difference. Within an orthonormal 
basis methodology, such as the R-matrix method we use here 
(Burke 2011; Berrington et al. 1995), the only way to account 
for the difference in orbitals is by introducing pseudo-orbitals 
(e.g., 3s, 3 p) such that the relaxed excited state can be described 
in terms of the ground state and the pseudo-orbitals via 

\s2s 2 2p A r = a ls2s g 2p g + c 2 ls2s 2 g 2ppp 

+ Cyls2s g 2p g 3s + ... . (4) 
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Figure 2. A 1 - matrix photoabsorption cross sections computed with and without relaxation effects and/or pseudo-resonance elimination. The red curve shows the 
definitive R-matrix calculation, which includes relaxation effects via the use of pseudo-orbitals. The blue curve in the upper plot shows results without relaxation, 
indicating a gross overestimate of the threshold energy position and cross section. The blue curve in the middle plot shows the results when using pseudo-orbitals but 
not using the pseudo-resonance elimination method (Gorczyca et al. 1995), giving large, unphysical features that permeate the threshold region and below. The two 
green curves show the IP asymptote and that reduced by 80% due to relaxation effects. The blue curve in the lower plot shows the results when including the additional 
ls2s 2 2p 3 3p and I s2s2p 4 3s target pseudo-states to give an approximate representation of the photoionization-excitation and double-photoionization channels, while 
eliminating all additional pseudo-resonances that are not associated with these pseudo-channels. 

(A color version of this figure is available in the online journal.) 


This procedure takes care of the relaxation effect and, in the 
present case, the reduction factor can be independently com- 
puted from simple multi-configuration Hartree-Fock (MCHF) 
calculations (Froese Fischer 1991) of |ci| 2 = 0.80. There is, 
thus, an analytically predicted reduction effect by a factor of 
0.80 due to relaxation, and the remaining 20% of the oscil- 
lator strength goes into photoionization-excitation and double 
photoionization. A discussion of these various contributions is 
given in the early experimental study of Ne i photoionization by 
Wuilleumier & Krause (1974; see, especially, their Figure 8). 

The effect of relaxation can also be seen by comparing the 
present //-matrix cross section with one where relaxation is not 
taken into account (Figure 2). The present, final //-matrix cross 
section asymptotically approaches the IP fit cross section of 
Verner et al. (1996), which is the scaled Is cross section after 
the latter has had the Is contribution multiplied by a factor 
of 0.80. This is the same overlap factor we compute from an 
independent MCHF calculation, thus independently confirming 
the 20% reduction effect. The original IP calculations (Reilman 
& Manson 1979), upon which Verner et al. (1996) based 
their fit, did not include relaxation effects. Consequently, their 
asymptotic value reflects the total photoabsorption cross section; 


this total includes shake-up and shake-off processes in addition 
to the direct Is photoionization (without secondary excitation 
or ionization). 

However, as also seen in Figure 2, even though the unrelaxed 
orbital results approach the full IP cross section (Verner et al. 
1996) asymptotically, they grossly overestimate the correct cross 
section just above the K-shell threshold; here, only the (relaxed) 
direct photoionization is energetically allowed. This overvalue 
is carried below threshold, leading to an unphysically enhanced 
resonance-oscillator strength. Moreover, the threshold energy 
position is also overestimated — by more than 10 eV — because 
of the inaccurate representation of the ls2s 2 2p 4 inner-shell 
vacancy state. 

It is, therefore, critical to account for relaxation effects; 
this is accomplished in the present theoretical methodology 
by including additional pseudo-orbitals in the atomic orbital 
basis set. However, without proper care, this procedure can 
lead to spurious pseudo-resonance structure (Gorczyca et al. 
1995), as is also shown in Figure 2. The present //-matrix re- 
sults, which were computed with pseudo-orbitals and with a 
proper elimination of pseudo-resonances (Gorczyca et al. 1995), 
are compared with similar //-matrix results without such an 
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Figure 3. Near- and above-threshold cross section: the red curve and black curve 
show the definitive /^-matrix and RMPS calculations, respectively. The two green 
curves show the IP asymptote and that reduced by 80% due to relaxation effects. 
The final fit formula is shown as the blue curve; the experimental results of Stolte 
et al. (1997), shifted by +0.58 eV, are given as the cyan data points; and the 
solid-state results of Henke et al. (1993) are given as the magenta squares. 

(A color version of this figure is available in the online journal.) 


devoid of pseudo-resonance structure; therefore, we choose this 
continuous data as the best representation of the cross section 
for energies above threshold. The R-matrix and RMPS cross 
sections just above threshold coincide with the data by Henke 
et al. (1993). Lastly, the measurements of Stolte et al. (1997) are 
consistent with the A 1 -matrix results throughout, as we address 
in Section 2.3.6. 


2.3.2. Analytical Fit to the 1j — »■ np Resonance Region 

The formula to fit the single-resonance photoabsorption cross 
section — parameterized by an absorption oscillator strength/, a 
resonance position E r , and a width T — is given by (see Bethe & 
Salpeter 1957, Equation (71.19)) 


OPA (E) = 


jt(k e e 2 )h df 
me dE ’ 


(5) 


where the oscillator strength per unit energy for an isolated 
resonance takes the form 


# = r/27r 

dE ' (E - E r ) 2 + ( r / 2) 2 ’ 


elimination. It is seen that the latter cross section exhibits large, 
spurious pseudo-resonance structure at higher energies. Fur- 
ther, these broad, unphysical resonance features permeate even 
down to the threshold region, resulting in an overestimate of 
the near-threshold cross section (and the resonance absorption 
oscillator strengths below threshold). By applying the pseudo- 
resonance elimination method (Gorczyca et al. 1995), the cross 
section becomes smooth throughout, and provides the most 
reliable resonance oscillator strengths as discussed in 
Section 2.3.1. The earlier A'-matrix calculations did not use 
a pseudo-resonance elimination method and, therefore, over- 
estimated the resonance absorption oscillator strengths (see 
Gorczyca & McLaughlin (2000) and their Figure 2). 

We can improve the asymptotic situation somewhat by includ- 
ing the orthogonal compliments to the \s2s 2 2p 4 0 + (relaxed) 
target states, namely the additional pseudo-states that are com- 
posed of the ls2s 2 2p 2 3p and I s2s2p A 3s configurations (with 
smaller »20% mixing of the \s2s 2 2p A configurations). This 
A 1 - matrix with pseudo-states (RMPS) method (Burke 2011), 
as implemented in the present codes following the develop- 
ments of Gorczyca & Badnell (1997), gives a somewhat crude, 
approximate description of the photoionization-excitation and 
double-photoionization channels, and, importantly, leads to the 
correct high-energy photoionization asymptote, as seen in the 
lower panel of Figure 2. The implementation of A'-matrix meth- 
ods on modern, massively parallel machines (Ballance & Griffin 
2006) will allow for a much larger, converged RMPS treatment 
of the problem. 

The findings thus far regarding the above-threshold cross 
section are summarized in Figure 3: Asymptotically, the present 
A'-matrix cross section approaches the IP results after the 
(dominant) Is contribution has been scaled by a factor of 0.8 
to account for relaxation, whereas the RMPS values show the 
correct asymptote but are still plagued by pseudo-resonances. 
The fit of Verner et al. (1996) at threshold is an extrapolation of 
the high-energy IP cross section of Reilman & Manson (1979), 
and, therefore, does not include the correct threshold rise as is 
seen in the A'-matrix results. On the other hand, the data of Henke 
et al. (1993), which are based on solid-state measurements, 
show the correct threshold and asymptotic cross sections and are 


i.e., it is equal to the discrete oscillator strength / times an 
energy-normalized Lorentzian, 

f r I2n 

J dE {E-E r ) 2 + (T/2) 2 ~ ' <?) 

Therefore, the photoabsorption profile can be characterized as 


&pa(E) = (if — 

(E 


T/2n 

F ,.) 2 + ( r / 2 ) 2 ’ 


( 8 ) 


with 


jtk e e 2 h 

P = — = 109.7626 Mb eV. 

me 

For an entire Rydberg series — characterized for each member 
by a principal quantum number «, resonance positions E n , 
width r (n-independent for inner-shell spectator Auger decay), 
and the oscillator strengths f„ — can be parameterized by a 
quantum defect ji, a threshold energy /*’ th , and an n-independent 
“strength”/): 


E 


n 


fn 


j-, ^ ^au 

th “ 2(« - /x) 2 ’ 

fo 

( n - M ) 3 ’ 


(9) 


where £ au = 27.21 leV. This discrete expression carries over 
to an analytic above-threshold cross section 


lim ctp\(E) = 



( 10 ) 


which must be considered when developing a consistent, con- 
tinuous formulation through threshold. 

Equation (9), based on quantum defect theoretical considera- 
tions, is precise in the limit n —*■ oo, but deficient for the lower 
resonances (n = 2, 3). The lowest members are more appropri- 
ately modeled by using separate (energy-dependent) quantum 
defects and oscillator strengths. 
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For multiple Rydberg series, if the interaction between them 
is neglected, the contribution from each series can be consid- 
ered separately. We apply this approach to the two dominant 
photoabsorption series in oxygen, namely \s2s 2 2p 4 { 4 P)np and 
\s2s 2 2p 4 ( L P)np (labeled, respectively, by the indices i s = 1 
and i s = 2), giving a two-series resonance cross section param- 
eterized as 


< S (E) = 


E 


'■=1 


«=«min=7 + l 


f0,n 

(n - p n ) 2 

rj2n 


(E - (El + Z 2 E m /(n - p n ) 2 )) + (r'72y 


fk 


0,oo 


Z 2 E ai 


' 1 1 

arctan 

2 IT 


El-E 

r/2 , 


(ii) 


The last term ensures that, since the below-threshold contri- 
bution has effectively been Auger broadened, i.e., convoluted 
with a Lorentzian of width T within each resonance energy 
interval A E ~ E au /(n — p) 2 , the step function, due to the 
above-threshold continuum photoionization, is likewise convo- 
luted near threshold: 
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where 0(E — £*) denotes the Heaviside step function at thresh- 
old. As the energy is increased above threshold, this expression 
is continuously extended to have the correct asymptotic tail 
as determined from our fit to the Henke et al. (1993) data in 
Equation (3), taking instead the form above threshold 
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In our R-matrix calculations, the oscillator strength (cross 
section) is found to be partitioned into the two dominant series 
by the fractions of 3/5 and 2/5 for the l.s2,y 2 2/? 4 ( 4 R)n/7 and 
\s2s 2 2p 4 ( 2 P)np series, respectively, instead of the statistical 
weighting of 2/3 and 1/3, because of channel coupling in 
the threshold region. As a result, the net oscillator strength 
density above threshold, which we find to be /o i00 = 0.132, is 
partitioned as = 0.6 hoo and = 0.4/ O ,oo- 

Our strategy is, first, to fit this expression to our present 
R-matrix results to get a good representation of the oscillator 
strengths and quantum defects, using the threshold energies and 
widths from the R-matrix runs. Then, the threshold energies are 
slightly adjusted and the analytical (Lorentzian) fit is further 
convoluted with the experimental (Gaussian) width to obtain a 
good fit to the experimental resonance spectrum. Before fitting 
this expression to our current results, we address resonance 
energy positions. 


2.3.3. Resonance Energy Positions from Astronomical Observations 

Oxygen K-shell photoabsorption in the ISM has been ob- 
served with both the Chandra and XMM-Newton satellite-borne 


Table 2 

Chandra Observations Used in This Work 


Source 

ObsID 

Date 

Exposure 

(ks) 

Read Mode 

4U 1636-53 

105 

1999 Oct 20 

29 

TE 


1939 

2001 Mar 28 

27 

TE 


6636 

2007 Jul 2 

26 

CC 


6635 

2006 Mar 22 

23 

CC 

4U 1735-44 

704 

2006 Jun 9 

24 

TE 


6637 

2006 Aug 17 

25 

CC 


6638 

2007 Mar 15 

23 

CC 

4U 1820-30 

1021 

2001 Jul 21 

9.6 

TE 


1022 

2001 Sep 12 

11 

TE 


6633 

2006 Aug 12 

25 

CC 


7032 

2006 Nov 5 

47 

CC 


6634 

2010 Oct 20 

26 

CC 

Cygnus X-l 

3407 

2001 Oct 28 

17 

CC 


3724 

2002 Jul 30 

8.8 

CC 

Cygnus X-2 

1102 

2003 Sep 23 

28 

TE 


8599 

2007 Aug 23 

59 

CC 


8170 

2007 Aug 25 

65 

CC 


10881 

2009 May 12 

66 

CC 

GX 9+9 

703 

2000 Aug 22 

20 

TE 


11072 

2010 Jul 13 

95 

TE 

XTEJ1817— 330 

6615 

2006 Feb 13 

18 

CC 


6616 

2006 Feb 24 

29 

CC 


6617 

2006 Mar 15 

47 

CC 


6618 

2006 May 22 

51 

CC 


Note. TE: TIMED; CC: CONTINUOUS. 


observatories. The High Energy Transmission Grating Spec- 
trometer (HETGS) of Chandra, in combination with the Ad- 
vanced CCD Imaging Spectrometer (ACIS), may provide the 
best spectral resolution with adequate sensitivity. It is exem- 
plified by the Juett et al. (2004) study, using the medium en- 
ergy gratings (resolving power of 0.023 A FWHM and an abso- 
lute wavelength accuracy of 0.01 1 A; Canizares et al. 2005) the 
HETGS to observe six Galactic X-ray sources. For the present 
work, we carry out a reanalysis of these observations in an at- 
tempt to improve the positions of the Oi Ka (l.y -* 2 p) and 
K/3 (l.y — 3 p) resonances. The XTE J1817— 330 source, pre- 
viously treated by Gatuzz et al. (2013), is included as well as 
some additional spectra for the sources considered by Juett et al. 
(2004). 

Observational specifications for the seven low-mass X-ray 
binaries used in this analysis are listed in Table 2. The ob- 
servations were taken in continuous clocking (CC) mode or 
time exposure (TE) mode. In CC mode, the temporal resolu- 
tion is increased to minimize the pileup effect (Cackett et al. 
2008). In TE mode, the ACIS instrument periodically reads the 
collected photons. All the spectrum files, response files, auxil- 
iary response files, and background files were taken from the 
Chandra Grating-Data Archive and Catalog TGCat. 1 1 We use 
the isis 12 package (version 1.6.2-18) for spectral fitting. 

We fitted all the observations for each source simultaneously 
using a simple power law+gaussians model in the oxygen- 
edge region (2 1-24 A). The power-law parameters were taken 
as independent, free parameters for each observation. We apply 
cash statistics because of the low signal-to-noise ratio in these 
spectra, which requires a minimal grouping of the spectra of 
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Table 3 

O i Kcr and K/i Line Positions (A) 


Source 

Key (Is 

-2 p) 



K/3 (Is 

- 3 p) 


Present 

Juett et 

al. 

(2004) 

Present 

Juett et al. (2004) 

4U 1636-53 

ii S0Q +0006 

23.507 

± 

0.011 

22.889f° 0 0 0 ° 0 5 4 

22.915 ±0.013 

4U 1735-44 

23.507 ± 0.009 

23.503 

± 

0.009 

22-890^«» 

22.861 ±0.006 

4U 1820-30 

23.509 ± 0.004 

23.514 

± 

0.010 

22.880!°d 0 0 0 i 9 o 

22.867 ± 0.015 

Cygnus X- 1 

23 507 +0003 

z,j.dv /_ 0004 

23.511 

± 

0.007 

oo qco+O.OII 
Zz.ooZ _ 0 0Q9 

22.8881° 0 ° 0 2 1 3 6 

Cygnus X-2 

23.508 ± 0.002 

23.508 

± 

0.004 

ZZ.ooj _ 0 006 

22.877™“ 

GX 9+9 

23 'SO3 +0006 

23.517 

± 

0.009 

22 QOO +0 008 
zz.yuu _ 0010 

22.906 ±0.018 

Mean position 3 

23.507 ± 0.005 

23.510 

± 

0.008 

22.886 ± 0.008 

22.8857 0 o°o 1 1 5 

XTEJ1817— 330 

23.506 ±0.001 




22.889 ± 0.004 


Mean position* 3 

23.507 ± 0.004 

23.510 

± 

0.008 

77 887+0.008 
ZZ.oo /_o 007 

22.8853°”, 


Notes. 

a Excluding XTE J1817— 330. 
b Including XTE J1817— 330. 


at least one count per spectral bin (see Humphrey et al. 2009; 
Baldi et al. 2012). Table 3 shows the Oi Ka and K/3 absorption 
line positions obtained from these fits. First, we list the mean 
position values excluding XTE J 1 8 17 — 330 to compare with 
those originally reported by Juett et al. (2004). When this source 
is taken into account, there is a slight decrease of the mean- 
value error bars. In the case of Cygnus X-l, the values in 
Juett et al. (2004) correspond to the average of ObsID 3407 
and ObsID 3742. From the confidence intervals, it may be 
appreciated that the present results are in agreement with the 
line positions estimated by Juett et al. (2004), but with improved 
statistics. 

The transition energy of the Oi Ka ( 1 ,v -» 2 p) line is also 
estimated using the well-exposed XMM-Newton RGS spectrum 
of Mrk 421 (Kaastra et al. 2006). This spectrum shows strong 
interstellar Ka absorption lines from O i at 23.5 138 ± 0.0022 A 
and Ovu at 21.6027 ± 0.0021 A. Since accurate theoretical 
values (Drake 1988; Cann & Thakkar 1992) and high-precision 
laboratory measurements (Engstrom & Litzen 1995) have been 
reported for the latter at 21.6015 A and 21.60195 ± 0.0003 A, 
respectively, we assume a small offset on the wavelength scale 
for this data set of 0.8 mA, which is well within the systematic 
uncertainty of the RGS. Correcting for this small difference 
(and assuming that the O vn line shows no intrinsic redshift), 
we find an energy of 527.30 ± 0.05 eV for the Oi Ka line. In 
summary, a list of the Oi Is — »■ 2 p and Is — > 3p line energies 
deduced from astronomical observations is given in Table 4. 
The listed Chandra line energies, unlike the XMM-Newton , 
were neither re-scaled with the Ovu Ka line nor Doppler 
corrected for the motion of the Earth around the Sun. Regarding 
the former issue, Gatuzz et al. (2013) quotes a line energy of 
21.593 ± 0.002 A for the Chandra observation of Ovu Ka 
toward XTE J 1817— 330, which is 9 mA short of the laboratory 
standard. If the wavelength scale is adjusted accordingly, our 
Chandra Oi Ka position in Table 4 would be reduced to 
527.26 ± 0.09 eV, and would be in better agreement with the 
XMM-Newton value. 

Recently, Liao et al. (2013) conducted a new and independent 
investigation of 36 Chandra HETG observations of 1 1 low-mass 
X-ray binaries, which accounts for the Galactic rotation velocity 
relative to the rest frame and uses a similar merging of corrected 
spectra. They use a Bayesian analysis to quantify systematic 
uncertainties and bias corrections, obtaining a resonance posi- 
tion with improved statistics. The resulting energy position of 


527.39 ± 0.02 eV, as listed in Table 4, is in agreement with our 
average observed value of 527.37 eV. 

2.3.4. Resonance Energy Positions from Laboratory Measurements 

We now consider the laboratory data for atomic oxygen. 
In Figure 4, we show the differences between the measured 
resonance energies for two experiments, namely Menzel et al. 
(1996) and Stolte et al. (1997). First, there is a systematic, 
almost linear change of the energy differences. This may be 
attributed to small remaining calibration uncertainties in at 
least one of the two datasets. Further, from the scatter, it is 
seen that the correlation between both datasets is much better 
than suggested by the formal error bars. This is probably 
because the error bars include a systematic uncertainty that 
may be nearly the same for all transitions. From a linear 
regression, we obtain for this energy difference (in eV units): 
A E = (36.357 ± 0.010) - (0.0670 ± 0.0016)£, with a scatter 
of 0.03 eV (much smaller than the nominal uncertainties of 
0.10 eV). Next, we compare the measurements of Krause (1994) 
and Caldwell et al. (1994) with those of Stolte et al. (1997). In 
this case, we find no significant slope (best fit 0.0029 ± 0.0037) 
and only a constant offset of 0.444 ± 0.028 eV. We conclude 
that the relative energy scales of Krause (1994), Caldwell et al. 
(1994), and Stolte et al. (1997) agree, and that, most likely, 
the energy scale of Menzel et al. (1996) is slightly off. All 
these datasets, however, show a different offset for their absolute 
energy scale. 

2.3.5. Resonance Energy Positions from Large MCHF Calculations 

The A- matrix calculations seek to span an indenumerably infi- 
nite number of bound, autoionizing, and continuum states of O i 
within a single, orthonormal basis of configurations and orbitals. 
It therefore becomes difficult to describe any specific state to a 
high degree of accuracy. In the present calculation, we are lim- 
ited to an active space of up to n = 2 physical orbitals and n = 3 
pseudo-orbitals. However, because the dominant Is -» 2p tran- 
sition energy is the source of a rather large (~ 0.5 -0.6 eV) 
discrepancy between observations and laboratory experiments, 
we can shed further light on the issue by appealing to sepa- 
rate, highly correlated theoretical calculations for the initial and 
final states. To this end, we use the sophisticated MCHF atomic 
structure package (Froese Fischer 1991) to perform a series of 
calculations using separate, large configuration-interaction (Cl) 
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Table 4 

Line Energies (eV) for O I 


Data Set 

Is -> 2 p 

Is -> 3 p 

A E 

Astronomical observations 
Chandra, average of seven sources 
XMM -Newton, Mrk 421 
Juett et al. (2004), six sources 

Average 

Chandra, Liao et al. (2013) 

527.44 ± 0.09 
527.30 ± 0.05 
527.41 +0.18 
527.37 

527.39 + 0.02 

541.72 + 0.18 
541.95 + 0.28 
541.77 + 0.40 

14.28 + 0.21 
14.65 + 0.33 
14.36 + 0.58 

Laboratory measurements 
McLaughlin et al. (2013) 

Stolte etal. (1997) 

Krause (1994), and Caldwell et al. (1994) 
Menzel et al. (1996) 

526.79 ± 0.04 
526.79 ± 0.04 
527.20 ± 0.30 
527.85 + 0.10 

541.19 + 0.04 

541.20 + 0.04 

541.27 + 0.15 

14.40 + 0.08 

14.41 ±0.08 

13.41 ±0.25 

MCHF calculations 

(«max = 6) 

527.49 




Note. The average of the Chandra and XMM data is found to be 527. 37 eV and is indicated in bold. 


Table 5 

MCHF Data for ls 2 2s 2 2p 4 ( 3 /n ls2s 2 2p 5 ( 3 P) in Oi 


^max 

E, 

(a.u.) 

Ef 

(a.u.) 

A E 
(eV) 

A 

fv 

2 

-74.85830 

-55.44337 

528.29 

0.133 

0.121 

3 

-74.99720 

-55.63645 

526.82 

0.107 

0.102 

4 

-75.06477 

-55.68599 

527.31 

0.098 

0.101 

5 

-75.08774 

-55.70510 

527.41 

0.093 

0.097 

6 

-75.09707 

-55.71152 

527.49 

0.097 

0.096 


Notes. Results are given as a function of « max , the maximum principal quantum 
number included in the active space expansion of configurations obtained by 
single and double promotions out of the initial or final configuration. Separate 
orbital bases are used for initial and final states, and relativistic corrections 
account for an additional ~0.03 eV to the transition energy, fa and fy are, 
respectively, the oscillator strengths in the length and velocity gauges. 

expansions, thus increasing the basis size to study the 
convergence of transition energies and oscillator strengths. 
Specifically, starting with the initial 1 j 2 2s 2 2/? 4 ( 3 P) configu- 
ration, we use a basis consisting of all configurations obtain- 
able by any single or double promotions from the outer n = 2 
orbitals into the active set of orbitals. For n max = 2, the ad- 
ditional ls 2 2s2p 5 and ls 2 2p 6 configurations are taken into 
account; for « max = 3, configurations such as 1 .v 2 2.v 2 2 /? 3 3 f 2 , 
ls 2 2s2p 4 3l, \s 2 2s 2 2p 2 3l3l' , \s 2 2s2p 2 3l3l' , and \s 2 2p 2 3l3l' 
(31 = {3s, 3p, 3d}) are also included. This procedure is re- 
peated for n max = 4,5,6, and, at each stage, a full-scale 
MCHF calculation is performed, optimizing each of the sep- 
arate orbitals from n = 1 to n = « max to produce a lengthy 
multi-configuration wavefunction for the initial state. This same 
procedure is repeated for the final state, re-optimizing all of 
the orbitals separately. Last, for given initial and final wave- 
functions (using completely different, non-orthogonal orbital 
bases), the absolute energies, transition energies, and oscil- 
lator strengths are computed. The results, listed in Table 5, 
show that the transition energy oscillates significantly between 
the observed and experimental values at first, but converges 
to a value consistent with that determined from the X-ray 
observations. 

2.3.6. Final Resonance Fit 

By fitting the expression in Equation (11) to the /(-matrix 
results, we obtain the parameters that are listed in Table 1. 
However, our initial fit used, in addition to the same widths as 


determined in the /(-matrix run (0.1348 eV and 0.1235 eV for 
series i s = 1 and i s — 2, respectively), the theoretical threshold 
energy positions E^ h = 544.74 eV and E 2 h = 549.67 eV. As 
is seen at all levels in Figure 5, this prescription provides 
an excellent fit to the /(-matrix results. The fit formula was 
next compared with the experimental cross section of Stolte 
et al. (1997) (shifted by +0.58 eV to position the Is — » 2 p 
resonance at 527.37 eV). However, to align our fit with these 
shifted experimental results, which comprise our assessment 
of the most accurate resonance positions, we had to shift our 
theoretical threshold energies by —0.2 eV and —0.35 eV for 
the two series i s = 1 and i s = 2, respectively. This has the 
simple effect of shifting each resonance of each series by these 
amounts. Further, to obtain the most meaningful comparison, 
it was necessary to convolute the fitting expression with a 
Gaussian of 182 meV FWHM to simulate the experimental 
resolution. Finally, it was necessary to upscale the n = 2 
oscillator strength to match the more reliable MCHF value of 
0.097 (see Table 5); the /(-matrix n = 2 oscillator strength 
was scaled down by a factor of 0.80 from the n — > oo 
series limit; and we use, as the final fit, the slightly increased 
value /q 2 = 0.867/q 1 ^, which ensures an oscillator strength 

of/ = /o,2/(2-/4) 3 ’ = 0-097- 

With these final adjustments, the resulting fit reproduces the 
experimental results well, except for the n = 2 resonance 
strength. In particular, the quantum defects for each of the 
two series, as determined from the fit to the /(-matrix results, 
align satisfactorily with the experimental values indicating that, 
regarding energy determination, the main source of error lies in 
determining accurate threshold positions. Further, with respect 
to the experimental spectrum, it seems that, if there is any 
error, it must be an overall global offset that we assume here 
to be —0.58 eV. The fit formula of Equation (11) could have 
just been modified by replacing the unit Lorentzian profiles 
by unit Voigt profiles, but it is more straightforward to simply 
convolute the resulting cross section numerically. Nevertheless, 
this highlights the flexibility of the analytical fitting formula 
in that the particular resonance shape can be, if desired, easily 
accounted for analytically. 

Finally, the experiment has a noticeable signal due to molec- 
ular O 2 contamination in the beam as evidenced by the strong 
Is tx* resonance at 531 eV. Consequently, there is an ad- 
ditional signal in the experiment throughout the K-edge re- 
gion, and the experimental procedure does not provide the most 
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Energy (eV) (Stolte et al. 1997) 


Figure 4. Comparison between the experimental resonance positions of Menzel 
et al. (1996) and Stolte et al. (1997). 


accurate benchmark away from resonance. On resonance, the 
signal is so strong and predominantly atomic in nature that 
the small molecular admixture would not affect the oscillator 
strength as much. 

3. COMPARISON WITH OTHER MODELS 

The three datasets of X-ray absorption currently used in 
spectral modeling codes we compare with the present fit are: (1) 
the R-matrix cross sections of Gorczyca & McLaughlin (2000); 
(2) thexSTAR modeling code, which uses the R-matrix cross 
sections of Garcia et al. (2005) except for independent data for 
the lowest Is —> 2 p resonance; and (3) spex, which is based 
on calculations with hfr (Cowan 1981) for resonances and data 
by Verner et al. (1996) elsewhere. A comparison is given in 
Figure 6 between these various approaches and the present fit, 
where several points must be brought to light. First, the cross 
section of Gorczyca & McLaughlin (2000), at least with respect 
to resonance positions, is close to the present fit since it is 
based on similar R-matrix calculations. However, the resonance 
oscillator strengths and the above-threshold cross section are 
higher, which we explain in Section 2.3.1 as being due to 
pseudo-resonance contamination present in the earlier R-matrix 
calculations. Further, the earlier R-matrix results show a minor 
discontinuity at the low-energy tail of the ls2s 2 2p 4 ( 4 P)3p 
resonance ( E ~ 538 eV) because of the sudden turn off of 
spectator Auger damping, which, as discussed in Section 2.3.1, 
we alleviate in the present study. 

The R-matrix results in xstar by Garcia et al. (2005) are 
seen in Figure 6 to be even higher than the present fit, 
regarding both the resonance oscillator strengths and the above- 
threshold cross section, and this is believed to be due to an 
insufficient treatment of both Cl and relaxation effects. In 
Garcia et al. (2005), relaxation is partially accounted for by 
optimizing the orbitals on a weighted sum of closed ls-shell and 
I ,v -vacancy states of O ii. Last, the spex model has gaps in the 
total oscillator strength density since only a finite number of 
terms are included for each Rydberg series. The unconvoluted 
natural widths are underestimated (in the original data, at least) 
since the spectator Auger decay was not taken into account; 
the included participator Auger width approaches zero as n 
increases (r„ ~ n 3 ). Additionally, the above-threshold cross 
section is matched to the IP results of Verner et al. ( 1 996), which, 
as discussed in Section 2.3.1, underestimate the threshold value. 



Figure 5. Fitting comparison: .K-matrix results (red curve), analytical fit formula 
to theoretical results (green curve), experimental results of Stolte et al. (1997), 
shifted by +0.58 eV (cyan data points), analytical fit formula using adjusted 
threshold energies and convoluted with the experimental resolution of 1 82 meV 
(cyan curve), Henke et al. (1993) data (magenta squares), Is — ► 2 p resonance 
position of 527.37 eV determined from observation (black vertical line, see text). 
(A color version of this figure is available in the online journal.) 


We also show the newer R-matrix results of McLaughlin et al. 
(2013), which are seen to closely reproduce the earlier R-matrix 
results of Gorczyca & McLaughlin (2000), including the over- 
estimate of the above-threshold cross section, as discussed in 
Section 2.3.1. One additional shortcoming in the new R-matrix 
results is that, since spectator Auger decay is not implicitly ac- 
counted for in that formulation, the predicted natural widths 
are underestimated (and, indeed, scale unphysically as 1/n 3 as 
n oo). 
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527 527.2 527.4 527.6 527.8 528 528.2 


Photon Energy (eV) 



Figure 6. Comparison between various datasets: the present fit model (red 
curve), SPEX (blue curve), /^-matrix results of Garcia et al. (2005; cyan curve), 
/^-matrix results of Gorczyca & McLaughlin (2000; green curve), /^-matrix 
results of McLaughlin et al. (2013; black curve), Henke et al. (1993) data 
(magenta squares), and present Chandra position of the Is — >- 2p resonance at 
527.37 eV (magenta vertical line). The earlier /^-matrix results of Gorczyca & 
McLaughlin (2000) and the more recent /^-matrix results of McLaughlin et al. 
(2013) are indistinguishable, except for the lowest Is — >- 2p resonance in the 
middle figure. 

(A color version of this figure is available in the online journal.) 

4. SUMMARY OF FITTING FORMULA 

The final expression for the photoabsorption cross section 
consists of the sum of the cross sections in Equation (1), where 
(T 2 S , 2 p{E), er[® s (£') and a^ ect (E) are given by Equations (2), (11), 
and (13), respectively, and the required fitting parameters are 


listed in Table 1. This final expression has several desirable 
features: 

1. It is an analytical formula easily transportable between 
different platforms and modeling codes; the Fortran routine 
used to generate a numerical photoabsorption cross section 
for all energies involves only about 100 lines of code. 

2. The formulation contains adjustable fitting parameters to 
best represent: (1) the K-edge positions; (2) the n -» oo 
energy-independent quantum defects and oscillator 
strengths; and (3) the energy-dependent quantum defects 
and oscillator strengths for the lower two resonances. From 
this fit, all relevant atomic parameters can be read off; for 
instance, the strongest Is -» 2p oscillator strength can be 
computed from our fit as 

/ = /o, 2 /(2 - /4) 3 = 0.097, 

and the integrated resonance strength is therefore given 
by fif = 10.65 Mb eV. Further modifications to these 
parameters can be made if so desired. 

3. The energy spectrum is optimized on the resonance posi- 
tions determined from a combined experimental and obser- 
vational assessment. 

4. A constant-resonance-width cross section — a Lorentzian 
profile that is predicted on physical grounds due to spec- 
tator Auger broadening — is implicitly included in the final 
expression, and can be further modified analytically to in- 
clude additional broadening effects. 

5. A consistent threshold formulation is obtained in that the 
lim^oo /„ for the (scaled) oscillator strength joins ana- 
lytically and smoothly with the above-threshold oscillator 
strength density df/dE. 

6. The consistent above-threshold cross section has the impor- 
tant factor of 0.80 reduction due to relaxation effects, and 
is then extended to higher energies to include the shake-up 
and shake-off processes. This results in photoionization- 
excitation and double-photoionization contributions to 
account for the 20% difference, giving the correct 
E — > oo high-energy asymptote, i.e., an accurate 
“shoulder.” 

5. CONCLUSION 

We have developed an analytical expression that encapsu- 
lates all of the important physics in X-ray absorption of atomic 
oxygen at all photon energies relevant to spectral modeling. 
For energies below or above the K-edge resonance region, we 
use simple parametric fits to our best assessment of the cross 
section based on a convergence between experimental and the- 
oretical data. The strong 1 .v — ► np resonances belonging to the 
two dipole-favored Rydberg series, on the other hand, require 
special attention regarding the oscillator strengths (and analytic 
continuation to the above-threshold direct Is cross section) and 
resonance positions. For this important region, we appeal to a 
combination of A 1 -matrix and MCHF theoretical calculations, 
laboratory experiments, and X-ray astronomical observations. 
An outstanding issue is the large discrepancy of ~0.6 eV be- 
tween several recent observational assessments and the latest 
laboratory experiments. Unconventionally, we choose to use the 
final calibration as suggested by the observations, since several 
sources and an independent large MCHF calculation tend to add 
credibility to this choice. Further, the recent laboratory exper- 
iments (Stolte et al. 1997; McLaughlin et al. 2013) calibrated 
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the photon energy scale using the molecular oxygen Rydberg 
resonance features, and it is unclear to us how accurately those 
molecular positions are known, especially considering the un- 
certainties we find in the atomic resonance positions. A repeat of 
those experimental measurements, calibrated instead to the more 
well-known CO and CO 2 K-edge features, will be performed in 
the near future (W. C. Stolte 2013, private communication), 
which will shed more light on the existing discrepancy. 

The ultimate goal of the present work is to establish a defini- 
tive, transparent, and portable photoabsorption cross section 
that can be incorporated in the two spectral modeling codes, 
namely xstar andsPEX. The consistent use of this developed 
photoabsorption expression in both methods to address molecu- 
lar abundances in the ISM, and clarify the existing controversy 
regarding the atomic-molecular fractions (Garcia et al. 2011), 
will be the subject of a subsequent follow-up paper. 
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